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Abstract. Dense, continuous pulsar timing observations over a 24-hr period provide a 
method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 
20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz 
Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), 
and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar 
observations to detect or constrain GWs typically at nanohertz frequencies. In the case of 
the IPTA’s nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713-1-0747, GW 
limits in the intermediate frequency regime can be produced. The negligible change in dispersion 
measure during the observation minimizes red noise in the timing residuals, constraining any 
contributions from GWs due to individual sources. At 10“® Hz, the 95% upper limit on strain 
is 10“^^ for GW sources in the pulsar’s direction. 


1. Introduction 

The direct detection of gravitational waves (GWs) is a goal for which pulsar timing arrays (PTAs) 
are well-suited. PTAs are collections of millisecond pulsars (MSPs) distributed throughout the 
sky that are regularly monitored with large, typically single-dish radio telescopes. MSPs are 
recycled neutron stars (NSs) with low spin-down rates. The rotational regularity of MSPs means 
that their beamed radio emission functions as a stable clock. Deviations from regularity may 
signify a GW-induced distortion in the spacetime metric, changing the apparent distance between 
pulses and thus changing the observed spin frequency of a pulsar. PTA-oriented scientific 
collaborations around the world (the EPTA collaboration pQ, NANOGrav [2], and the PPTA 
collaboration [H] ; working together as the IPTA collaboration i) have been regularly monitoring 
their respective sets of pulsars for over a decade. The best-case GW detection would show the 
correlated signature between the pulse times-of-arrival (TOAs) from different pulsars predicted 
in [5]. Such a signal would be free of spurious effects along a particular line-of-sight (LoS), such 
as long-term, stochastic rotational fluctuations in a particular pulsar, also known as timing noise 
or spin noise. Although GW detection has not yet occurred with PTAs, constraints on galaxy 
evolution models have already been obtained iiEii, in some cases limiting the presence of GWs 
using multiple-pulsar data, and in others using data from a single pulsar. These non-detection 
limits are at frequencies near a nanohertz. 

2. Data Reduction, Pipeline, and Limits on Individnal Sources 

The PSR J1713-I-0747 24-hour global campaign [9] was a day-long, continuous observation using 
nine radio telescopes on 22-Jun 2013 (MJD 56465 - 56466). PSR J1713-I-0747 was chosen 
because of the low root-mean-squared (rms) timing residuals with respect to its timing model, 
and because it is regularly observed by all the PTA collaborations. The Arecibo, Effelsberg, 
GMRT, Green Bank, LOFAR, Lovell, Nangay, Parkes, and WSRT radio telescopes were all used 
for the observation. For the much longer datasets usually used in PTAs, red noise in dispersion 
measure (DM) variations propagates into timing residuals |T0]. DM is the quantity that measures 
the integrated column density of free electrons along the LoS to a pulsar. The 24-hr campaign 
obtained high signal-to-noise (S/N) TOAs at 1.4 GHz, showing that the timing residuals were 
free of any red noise due to DM variations. Noise in the residuals due to interstellar scintillation 
(ISS) was also modeled and shown to be white. This procedure was performed in [9j using 
the frequency-dependent (ED) parameters [H] that modeled the pulse profile’s evolution with 
frequency. Otherwise, performing a weighted sum will generate non-white noise structures due 
the the fact that ISS creates a varying S/N across the observing band. In that case, the varying 
S/N would result in different TOAs in each sub-band due to the intrinsic profile’s changing 
shape with frequency. 
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The timing residuals across the varying observation bandwidths are shown as a weighted 
sum, which we refer to here as the “broadband” residuals, in Fig. 1. Around 8 hr, ISS was 
occurring, which decreased the S/N. The lower S/N ratio was in spite of the fact that the 
GMRT’s sensitivity is comparable to that of other telescopes. These otherwise high-S/N data 
motivate a search for high-frequency GWs from individual (or continuous wave; GW) merging 
supermassive black hole binaries (SMBHBs), rather than a stochastic background of GWs from 
such sources from 10“® - 10“^ Hz. For the 24-hr campaign, we are sensitive to GWs at 10“® - 
10“^ Hz. At these high GW frequencies, we do not expect sensitivity to the GW background 
because of its powerlaw spectrum, which dominates at low frequencies. However, a nearby 
SMBHB with a sufficient chirp-mass could produce GWs. jT2] also demonstrates that the 10“^ 
" 10“^ Hz range may be a good discovery band for unanticipated, periodic GW sources. We use 
the same detection procedure for sinusoidal GWs as m, where details of the pipeline can be 
found. The only difference between the pipeline as typically used for a PTA dataset and for the 
present dataset is simply that only one pulsar, PSR J1713-I-0747, is included in the array. [T3] 
contains several detection algorithms; here we utilize the frequentist, fixed-noise model, given 
that we have already established the white-noise character of the dataset. The pipeline accounts 
for the uneven sampling and the varying S/N of the data over time. An additional pulse phase 
term due to any GW perturbations on the pulsar itself is included. 
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Figure 1. The 1.4 GHz broadband, 8- 
telescope timing residuals for the 24-hour 
global campaign on pulsar J1713-I-0747. 
The white-noise character of these residuals 
motives a search for GW in the frequency 
range 10“® - 10“^ Hz. The integration time 
is 120 s, the variation in rms is due to ISS 
and changing observation bandwidth, and the 
error bars are the 1 a residual uncertainty. 



Figure 2. GW limits from the J1713 24- 
hr global campaign (red) from individual GW 
sources in random sky directions. The dashed 
curve shows the GW limit in the direction 
of the pulsar. The other GW limits shown 
are: pulsar J0437 [Hj (blue), pulsar B1937 
[15] (magenta), NANOGrav [T3] (black), the 
EPTA [16] (green), and the PPTA |T7] (cyan). 


We obtain two 95% upper limit non-detection curves (Fig. 2) for the 24-hour global data: 
one for individual sources in a random direction in the sky (solid red curve), and the other for 
sources in the direction of PSR J1713-I-0747 itself (dashed red curve). Both show limits in GW 
strain h as a function of GW frequency (/gw)- We chose a 0.01 radian cone for the directional 
limit. A spin-down rate fit is also included. The directional limit curve, which has h ~ 10“^^ 
at 10“^ Hz, provides about a factor of 10 improvement. Also plotted are the single-pulsar high- 
frequency limit curves from [H] and [15], from PSR J0437-4715 (magenta) and PSR B1937-I-21 
(blue), respectively, for GW sources from a random sky direction. The random-direction GW 



























































































































limit curve from [T3] (black) is also shown, even though it is derived from 17 NANOGrav pulsars 
rather than from a single pulsar. The curve is included because it was generated from the same 
CW detection pipeline; many-pulsar CW limits from the EPTA [T^ (green) and the PPTA 
|17j (cyan), both random-direction, are also shown. We note that the NANOGrav curve can 
be extrapolated to meet the upper curve from the 24-hr campaign. This is consistent with 
the fact that PSR J1713-I-0747 influences the GW search procedure across all pulsars because 
of its exceptionally high residual S/N. Both |18j and the Cassini radio science investigation 
[19] also derived limits on GWs in similar /gw bands, but these apply only to a stochastic 
background and we do not consider them here. There are no known GW sources near this band 
that are strong enough to have been detected. According to |20|, the double pulsar J0737-3039 
1.1 kpc away [2Tl[22l[23] should produce hmax ~ 10“^^ at /gw ~0.2mHz. Indeed, the current 
GW frequency of J0737-3039 is the highest known amongst NS binaries, and is closest to the 
GW frequency band considered in this work. Further considerations about other potential GW 
sources will be explored in a future publication. Meanwhile, we recommend further long-track 
studies, especially those optimized for different sky directions, or those coordinated between 
different telescopes monitoring different pulsars simultaneously. 
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